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Abstract. An investigation of classical chaos and quantum chaos in gauge fields and fermion 
fields, respectively, is presented for (quantum) electrodynamics. We analyze the leading Lya- 
punov exponents of U(l) gauge field configurations on a 12^ lattice which are initialized by 
Monte Carlo simulations. We find that configurations in the strong coupling phase are sub- 
stantially more chaotic than in the deconfinement phase. Considering the quantum case, 
complete eigenvalue spectra of the Dirac operator in quenched Ad compact QED are studied 
on 8'' X 4 and 8^ X 6 lattices. We investigate the behavior of the nearest-neighbor spacing 
distribution P(s) as a measure of the fiuctuation properties of the eigenvalues in the strong 
coupling and the Coulomb phase. In both phases we find agreement with the Wigner surmise 
of the unitary ensemble of random-matrix theory indicating quantum chaos. 

1. Lyapunov exponents in Minkowskian U(l) gauge theory. 

1.1. Classical chaotic dynamics from Monte Carlo initial states. Cha- 
otic dynamics in general is characterized by the spectrum of Lyapunov expo- 
nents. These exponents, if they are positive, reflect an exponential divergence 
of initially adjacent conflgurations. In case of symmetries inherent in the 
Hamiltonian of the system there are corresponding zero values of these expo- 
nents. Finally negative exponents belong to irrelevant directions in the phase 
space: perturbation components in these directions die out exponentially. Pure 
gauge flelds on the lattice show a characteristic Lyapunov spectrum consisting 
of one third of each kind of exponents T]. Assuming this general structure of 
the Lyapunov spectrum we investigate presently its magnitude only, namely 
the maximal value of the Lyapunov exponent, Lmax- 

The general deflnition of the Lyapunov exponent is based on a distance 
measure d(t) in phase space, 

L := lim lim -In^^. (1) 

t^ood(o)^o t d{0) 

In case of conservative dynamics the sum of all Lyapunov exponents is zero 
according to Liouville's theorem, ^ Lj = 0. We utilize the gauge invariant 
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distance measure consisting of the local differences of energy densities between 
two 3d field configurations on the lattice: 

i^Ep Itrf/p-trf/;|. (2) 

Here the symbol Y2p stands for the sum over all Np plaquettes, so this distance 
is bound in the interval (0,2A^) for the group SU(N). Up and Up are the 
plaquette variables, constructed from the basic link variables U^^i, 

U^,i = exp {aAl^,r) , (3) 

located on lattice links pointing from the position x = {xi,X2,X3) to a; + acj. 
The generators of the group are T'^ = —igT'^/2 with t'^ being the Pauli matrices 
in case of SU(2) and A^^ is the vector potential. The elementary plaquette 
variable is constructed for a plaquette with a corner at x and lying in the ij- 
plane as Ux,ij = Ux,iUx+i,jUl_^_j^^Ulj. It is related to the magnetic field strength 

Ux,^j = exp {eijkaB^ f^T") . (4) 
The electric field strength ^ is related to the canonically conjugate momen- 
tum Px^i = Ux,i via 

= ^tr {t<^Ux,uI) ■ (5) 

The Hamiltonian of the lattice gauge field system can be casted into the 
form 

^■{p,p) + i-i(r/,v) 



. - (6) 

Here the scalar product stands for {A,B) = ^ti^AB"^). The staple variable 
is a sum of triple products of elementary link variables closing a plaquette 
with the chosen link U. This way the Hamiltonian is formally written as a 
sum over link contributions and V plays the role of the classical force acting 
on the link variable U. 

Initial conditions chosen randomly with a given average magnetic energy per 
plaquette have been investigated in past years We prepare the initial field 
configurations from a standard four dimensional Euclidean Monte Carlo pro- 
gram on a 12^ X 4 lattice varying the gauge coupling g [3] . We relate such four 
dimensional Euclidean lattice field configurations to Minkowskian momenta 
and fields for the three dimensional Hamiltonian simulation by identifying a 
fixed time slice of the four dimensional lattice. 

1.2. Chaos, confinement and continuum limit. We start the presenta- 
tion of our results with a characteristic example of the time evolution of the 
distance between initially adjacent configurations. An initial state prepared 
by a standard four dimensional Monte Carlo simulation is evolved according to 
the classical Hamiltonian dynamics in real time. Afterwards this initial state 
is rotated locally by group elements which are chosen randomly near to the 
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Figure 1. Exponentially diverging distance in real time of 
initially adjacent U(l) field configurations on a 12'^ lattice pre- 
pared at /? = 0.9 in the confinement (left) and at /3 = 1.1 in the 
Coulomb phase (right). 



unity. The time evolution of this slightly rotated configuration is then pur- 
sued and finally the distance between these two evolutions is calculated at the 
corresponding times. A typical exponential rise of this distance followed by a 
saturation can be inspected in Fig.^from an example of U(l) gauge theory for 
two values of /3 = 1/(7^ in the confinement phase and in the Coulomb phase. 
While the saturation is an artifact of the compact distance measure of the 
lattice, the exponential rise (the linear rise of the logarithm) can be used for 
the determination of the leading Lyapunov exponent. The left plot exhibits 
that in the confinement phase the field has larger Lyapunov exponents than 
in the Coulomb phase shown in the right plot. 

The main result of the present study is the dependence of the leading Lya- 
punov exponent I/max on the inverse coupling strength P, displayed in Fig. |21 
As expected the strong coupling phase is more chaotic. The transition reflects 
the critical coupling to the Coulomb phase. 

An interesting result concerning the continuum limit can be viewed from 
Fig. 121 which shows the energy dependence of the Lyapunov exponents for the 
U(l) theory. One observes an approximately quadratic relation in the weak 
coupling regime. From scaling arguments one expects a functional relationship 
between the Lyapunov exponent and the energy H] 

L{a) oc a''-^E''{a), (7) 

with the exponent k being crucial for the continuum limit of the classical field 
theory. A value of A; < 1 leads to a divergent Lyapunov exponent, while k > 1 
yields a vanishing L in the continuum. The case A; = 1 is special allowing for 
a finite non-zero Lyapunov exponent. Our analysis of the scaling relation ((7j) 
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Figure 2. Lyapunov exponents of 100 U(l) field configurations 
as a function of coupling. 
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Figure 3. Average maximal Lyapunov exponent as a function 
of the scaled average energy per plaquette ag^E. The U(l) gauge 
theory shows an approximately quadratic behavior in the weak 
couphng regime. 

gives evidence that the classical compact U(l) lattice gauge theory has k ^ 2 
and with L(a) — > a regular continuum theory. 

2. Quantum chaos in compact Euclidean QED. 

2.1. Manifestation of quantum chaos. The fluctuation properties of the 
eigenvalues of Dirac operator for quantum chromodynamics (QCD) on a lat- 
tice in Euclidean space-time have attracted much attention in the past few 
years. In Ref. |^ it was first shown for SU(2) lattice gauge theory that certain 
features of the spectrum of the Dirac operator are described by random-matrix 
theory (RMT). In particular the so-called nearest-neighbor spacing distribu- 
tion P{s), i.e. the distribution of the spacings s of adjacent eigenvalues on 
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the "unfolded" scale, agrees with the Wigner surmise of RMT. According to 
the Bohigas-Giannoni-Schmit conjecture |.61, quantum systems whose classical 
counterparts are chaotic have a P{s) given by RMT whereas systems whose 
classical counterparts are integrable obey a Poisson distribution P{s) = e~*. 
Therefore, the specific form of P{s) is often taken as a criterion for "quantum 
chaos". However, there is no accepted proof of the Bohigas-Giannoni-Schmit 
conjecture yet. The field of quantum chaos is still developing and there are 
many open conceptual problems I7J. Applying this conjecture it was recently 
demonstrated that QCD is chaotic, both in the confinement and the quark 
gluon plasma phase [H]. 

A number of interesting results have been established for chaotic dynam- 
ics in classical gauge theories. Lattice gauge theories are chaotic as classical 
Hamiltonian dynamical systems [2]. Furthermore, it was found that the lead- 
ing Lyapunov exponent of SU(2) Yang-Mills field configurations indicates that 
configurations corresponding to the deconfinement phase are chaotic although 
they are less chaotic than in the strong coupling phase at finite temperature 
P] . The scaling of the maximal Lyapunov exponent in the classical continuum 
limit was studied in Ref. It was suggested that Abelian gauge theories be- 
have regularly in the continuum limit whereas non- Abelian gauge theories are 
chaotic in the continuum, although the exact scaling relation is still an open 
problem. Chaos to order transitions were observed in a spatially homogeneous 
SU(2) Yang-Mills-Higgs system and in a spatially homogeneous SU(2) Yang- 
Mills Chern-Simons Higgs system jHlE]- In Ref. [H] a chaos to order transition 
was also seen on the quantum level, i.e. a smooth transition from a Wigner to 
a Poisson distribution was found. A transition in P{s) from Wigner to Poisson 
behavior was further observed at the metal-insulator transition of the Ander- 
son model JT]. Further, the suppression of the characteristic manifestations 
of dynamical chaos by quantum fluctuations was analyzed in the context of 
spatially homogeneous scalar electrodynamics and for a + 1-dimensional 
space-time A^-component 0^ theory in the presence of an external field |13]. 
These chaos to order transitions were seen in spatially homogeneous models 
and not for the full classical field theory. The relationship to properties of the 
quantum field theory is an interesting issue. 

Here we focus on the Dirac operator for quenched 4d compact quantum elec- 
trodynamics (QED) to search for the possible existence of a transition from 
chaotic to regular behavior in Abelian lattice gauge theories. In particular, we 
are interested in the nearest-neighbor spacing distribution of the eigenvalues 
of the Dirac operator across the phase transition from the strong coupling to 
the Coulomb phase. In the strong coupling region Abelian as well as non- 
Abelian lattice gauge theories are in a confined phase ;14J. For compact QED 
this means that for couplings /5 < /?c ~ 1-01 the electron is confined. How- 
ever, when crossing the phase transition the conventional Coulomb phase is 
observed. It is an interesting question if the difference between the Coulomb 
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phase in QED and the quark-gluon plasma phase in QCD has an influence on 
the level repulsion of the corresponding Dirac spectra. 

2.2. Quantum chaos of fermion fields. We generated gauge field config- 
urations using the standard Wilson plaquette action for U(l) gauge theory, 

5G(^/)=/35^^(l-cosep), (8) 

where Ui = U^^^ = exp(i02;,At), with Qx,fi £ [~7r,7r), are the field variables 
defined on the links / = (x, /i). The plaquette angles are 9p = 9^^^ + Qx+fi,u — 
Qx+u^i — Qx,u- We simulated 8'^ x 4 and 8^ x 6 lattices at various values of the 
inverse gauge coupling P = 1/ both in the strong coupling and the Coulomb 
phase. Typically we discarded the first 10000 sweeps for reaching equilibrium 
and produced 20 independent configurations separated by 1000 sweeps for 
each j3. Because of the spectral ergodicity property of RMT one can replace 
ensemble averages by spectral averages ^3] if one is only interested in the bulk 
properties. Thus a few independent configurations are sufficient to compute 
P{s). 

On the Euclidean lattice the Dirac operator Ip = p + ig for staggered 
fermions 

Mx,x' = 2 X^^^M {^x+fi,x'Ux,f, - Sx-[,,x'Ul^f) (9) 

is anti-Hermitian so that all eigenvalues are imaginary. For convenience we 
denote them by iXn and refer to the A„ as the eigenvalues in the following. 
Because of {p, 75} = the occur in pairs of opposite sign. All spectra were 
checked against the analytical sum rules 

5^A„ = and Y,\l = V, (10) 

n A„>0 

where V is the lattice volume. 

In RMT one has to distinguish between different universality classes which 
are determined by the symmetries of the system. So far the classification 
for the QED Dirac operator has not been done. Our calculations show that 
in the case of the staggered 4(i compact QED Dirac matrix the appropriate 
ensemble is the unitary ensemble. Although from a mathematical point of 
view this is the simplest one, the RMT result for the nearest-neighbor spacing 
distribution is still rather complicated. It can be expressed in terms of so- 
called prolate spheroidal functions, see Ref. [12] where P{s) has also been 
tabulated. A good approximation to P{s) is provided by the Wigner surmise 
for the unitary ensemble 

P(.) = ^.^e-^\ (11) 

We have simulated 8^ x 4 lattices at /5 = 0, 0.90, 0.95, 1.00, 1.05, 1.10, 1.50 
and 8=^ X 6 lattices at /3 = 0.90, 1.10, 1.50. All results are similar to those 
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Figure 4. Nearest-neighbor spacing distribution P{s) of the 
Dirac operator for compact U(l) theory in the strong cou- 
phng phase for P = 0.90 (left) and in the Coulomb phase for 
/3 = 1.10 (right). The histogram represents the lattice data on 
an 8^ X 6 lattice averaged over 20 independent configurations. 
The full curve is the Wigner distribution of Eq. (fTT|l for the uni- 
tary ensemble of RMT. For comparison the Poisson distribution 
P{s) = e~* is also indicated by the dashed line. 




Figure 5. Nearest-neighbor spacing distribution P{s) of the 
analytically calculated eigenvalues of Eq. ()12|) for a free Dirac 
operator on a 53 x 47 x 43 x 41 lattice (histogram) compared 
with the Poisson distribution P{s) = (solid line). 

selected for the plots. The left plot in Fig. |3] shows the nearest-neighbor spac- 
ing distribution P{s) for [3 = 0.90 in the confined phase averaged over 20 
independent configurations on the 8^ x 6 lattice compared with the Wigner 
surmise for the unitary ensemble of RMT of Eq. ()11|) . Good agreement is 
found. According to the Bohigas-Giannoni-Schmit conjecture this means the 
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system can be regarded as chaotic in the strong couphng region. The right plot 
in Fig. 0] shows the nearest-neighbor spacing distribution P{s) for (3 = 1.10 
in the Coulomb phase again averaged over 20 independent configurations and 
compared with the Wigner surmise ()11|). The agreement of the lattice data 
with the RMT predictions is interpreted gnal that quantum chaos sur- 

vives the phase transition. We find no deviation up to the maximum coupling 
considered, P = 1.50. 

In the strong coupling phase the result holds down to /3 = 0. Therefore, we 
tend to interpret our, as well as previous |H1E1, results in the sense that the dis- 
order of the gauge field configurations [21 E] is responsible for the chaotic char- 
acteristics of the spectrum of the Dirac operator. In contrast to that: The free 
fermion theory is non-chaotic and the corresponding nearest-neighbor spacing 
distribution obeys a Poisson distribution. This is illustrated in Fig. El where 
P{s) is obtained from the analytical eigenvalues of the free Dirac operator on 
a 53 X 47 X 43 X 41 lattice: 



Here a is the lattice constant, is the number of lattice sites in /i-direction, 
and = 0, — 1. We used an asymmetric lattice with being primes 
and restricted the range to (L^ — l)/2 instead of — 1 in each direction to 
avoid degeneracies of the free spectrum. 

3. Conclusion. In the underlying article we performed a comparison of clas- 
sical chaos and quantum chaos in fundamental field theories of physics, exem- 
plified for the U(l) theory of electrodynamics. This is not a direct comparison, 
however, since it deals with the gauge field in classical theory and with fermions 
in the quantum case. It turned out that the classical U(l) field is chaotic in the 
confinement phase with decreasing Lyapunov exponents towards the Coulomb 
phase. A scaling analysis indicates a regular continuum theory as one ex- 
pects from the Maxwell equations. On the other hand, our investigation of 
the quantized fermion field fulfills the criterion for quantum chaos both in the 
confinement and Coulomb phase. A scaling analysis was not possible for the 
quantum case (due to the lack of a /3-function) which could cover the transi- 
tion to a regular theory. Nevertheless, the free Dirac operator, in absence of a 
covariant derivative and minimal gauge coupling, exhibits regular behavior. 

It would be interesting to study the direct counterpart of the classical gauge 
field after quantization. A similarly accurate determination of the eigenvalue 
spectrum of the gauge sector necessitates to construct the corresponding Fock 
space and to diagonalize high- dimensional matrices which seems to be out of 
reach for 4d QED/QCD. On the other hand, chaos studies of the classical limit 
of the fermion field would also be of interest but have not yet been attempted. 
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